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SUMMARY 

The problem of controlling large, flexible space systems has been the 
subject of considerable research. Many approaches to control system 
synthesis have been evaluated using computer simulation. In several cases, 
ground experiments have also been used to validate system performance under 
more realistic conditions. There remains a need, however, to test 
additional control laws for flexible spacecraft and to directly compare 
competing design techniques. In this paper an NASA program is discussed 
which has been initiated to make direct comparisons of control laws for, 
first, a mathematical problem, then an experimental test article is being 
assembled under the cognizance of che Spacecraft Control Branch at the NASA 
Langley Research Center with the advice and counsel of the IEEE Subcom- 
mittee on Urge Space Structures. The physical apparatus will consist of a 
softly supported dynamic model of an antenna attached to the Shuttle by a 
flexible beam. The control objective will include the task of directing 
the line-of-sight of the Shuttle/antenna configuration toward a fixed 
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target, under conditions of noisy data, limited control authority and 

random disturbances. The open competition started in the early part 

of 1984. Interested researchers are provided information intended to 

facilitate the analysis and control synthesis tasks. A workshop is planned 

for early December at the NASA Langley Research Center to discuss and 
compare results. 


INTRODUCTION 


y future spacecraft will be large and consequently quite flexible. 
As the size of antennae is increased, the frequencies of the first flex- 
ible modes will decrease and overlap the pointing system bandwidth. It 
will no longer be possible to use low gain systems with simple notch 
filters to provide the required control performance. Multiple sensors and 
actuators, and sophisticated control laws will be necessary to ensure 

stability, reliability and the pointing accuracy required for large, 
flexible spacecraft. 

Control of such spacecraft has been studied with regard given to 
modeling, order reduction, fault management, stability and dynamic system 
performance. Numerous example applications have been used to demonstrate 
specific approaches to pertinent control problems. Both computer simula- 
tions and laboratory experiment results have been offered as evidence of 
the validity of the approaches to control large, flexible spacecraft. 
Concerns remain, however, because of the chronic difficulties in control- 
ling these lightly damped large-scale systems. Because of these concerns 
and because of the desire to offer a means of comparing technical 
approaches directly, an NASA/IEEE Design Challenge is being offered. An 
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experimenCal test article is being assembled under the cognizance of the 
Spacecraft Control Branch at the NASA Langley Research Center with the 
| advice and counsel of the IEEE (COLSS) Subcommittee on Urge Space 

Structures. This Spacecraft Control Laboratory Experiment (SCOLE) will 
t serve as the focus of a design challenge for the purpose of comparing 

directly different approaches to control synthesis, modeling, order 
reduction, state estimation and system identification. 

The configuration of the SCOLE will represent a large antenna attached 
to the Space Shuttle orbiter by a flexible beam. This configuration was 
chosen because of its similarity to proposed space flight experiments and 
proposed space-based antenna systems. This paper will discuss the "Design 
Challenge" in terms of both a mathematical problem and a physical experi- 
mental apparatus. The SCOLE program is not part of any flight program. 

SYMBOLS 

a acceleration vector ft/sec^ 

A beam cross section area 

c observation matrix 

d noise contaminating direction cosine matrix measurements 

e line-of-sight error 

E modulus of elasticity 

f concentrated force expressions 

F4 force vector 

g concentrated moment expressions 

„ GI torsional rigidity 

I moment of inertia matrix for entire Shuttle/antenna configuration 
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moment of inertia matrix. Shuttle body 

moment of inertia matrix, reflector body 

beam cross section moment of inertia, roll bending 

beam cross section moment of inertia, pitch bending 

beam polar moment of inertia, yaw torsion 

length of the reflector mast, beam 

control moment applied to the Shuttle body 

control moment applied to the reflector body 

disturbance moment applied to the Shuttle body 

mass of entire Shuttle/antenna configuration 

mass of Shuttle body 

mass of reflector body 

mass density of beam 

beam position variable 

direction cosine matrix, Shuttle body () = T () 

earth 1 Shuttle body 

direction cosine matrix, reflector body () = T () 

earth 4 reflector 

inertial velocity. Shuttle body b ° dy 

inertial velocity, reflector body 

lateral deflection of beam bending in y-z plane 

lateral deflection of beam bending in x-z plane 

angular deflection of beam twisting about z axis 

position variables 

displacement of proof-mass actuator 

line-of-sight pointing requirement 

noise contaminating angular velocity measurements 


t 


1 
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0,<t»,y pitch, roll, heading 

C damping ratio 

4 t noise contaminating acceleration measurements 

“l angular velocity of Shuttle body 

t o >4 angular velocity of reflector body 

DISCUSSION 

The objective of the NASA- IEEE Design Challenge concerning the control 
of flexible spacecraft is to promote direct comparison of different 
approaches to control, state estimation and systems identification. The 
design challenge has principal parts, the first using a mathematical model, 
and the second using laboratory experimental apparatus. The specific parts 

of the Spacecraft Control Laboratory Experiment (SCOLE) program will be 
discussed in detail. 

Control Objectives 

The primary control task is to rapidly slew or change the line-of- 
sight of an antenna attached to the space Shuttle orbiter, and to settle or 
damp the structural vibrations to the degree required for precise pointing 
of the antenna. The objective will be to minimize the time required to 
slew and settle, until the antenna line-of-sight remains within the 
angle 6. A secondary control task is to change direction during the 
"on- target" phase to prepare for the next slew maneuver. The objective is 
to change attitude and stabilize as quickly as possible, while keeping the 
line-of-sight error less than 6. 
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Math Model Dynamics 


The initial phase of the design challenge will use a mathematical 
model of the Shuttle orbiter/antenna configuration. It is necessary to 
obtain a balance, of course, between complex formulations which might be 

more accurate and simplified formulations which ease the burden of 
analysis. 

The dynamics are described by a distributed parameter beam equation 
with rigid bodies, each having mass and inertia at either end. One body 
represents Space Shuttle orbiter; the other body is the antenna reflector. 
The equations for the structural dynamics and Shuttle motion are formed by 
adding to the rigid-body equations of motion, beam-bending and torsion 
equations. The boundary conditions at the ends of the beam contain the 
forces and moments of the rigid Shuttle and reflector bodies. The 
nonlinear kinetmatics couples the otherwise uncoupled beam equations. 
Additional terms represent the action of two, 2-axis proof-mass actuators 
at locations on the beam chosen by the designer. 

The rigid-body equations of motion for the Shuttle body are given by: 

“1 ' - ^‘‘vri + ", * % * y,> 
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Similarly, for the reflector body. 
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The direction cosine matrices defining the attitudes of the Shuttle and 
reflector bodies are given by: 


•T ~ T 
T = - u) T 
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*T ~ T 
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The direction cosine matrices defining the attitudes of the Shuttle and the 
reflector bodies are related to the beam end conditions. 
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The equations of motion for the flexible beam-like truss connecting the 
reflector and Shuttle bodies consist of standard beam bending and torsion 
partial differential equations with energy dissapative terms which enable 
damped modes with constant characteristics for fixed, though dynamic, end 
conditions. The system of equations can be viewed as driven by changing 

end conditions and forces applied at the locations of the proof-mass 
actuators. 

ROLL BEAM BENDING: 
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PITCH BEAM BENDING : 
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YAW BEAM TORSION: 
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{PROOP-MASS ACTUATOR} 
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{REFLECTOR BODY FORCE} 
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{SHUTTLE BODY, MOMENTS} 


g 0,2 = 0 


(PROOF-MASS ACTUATOR, MOMENT} 


g 0,3 = 0 


(PROOF-MASS ACTUATOR, MOMENT} 
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{reflector body, moment} 


The angular velocity of the reflector body is related to the Shuttle 


O 13 c O 
R b = -130 o o 
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The line-of-sight error described in figure 2 is affected by both the 
pointing error of the Shuttle body and the misalignment of the reflector 
due to the deflection of the beam supporting the reflector. The line-of- 
sight is defined by a ray from the feed which is reflected at the center of 
the reflector. Its direction in the Shuttle body coordinates is given by: 


where : 


LOS 



[ R I<\ - v • ■>*] 

IK - “p- 2 

r a r r ~ r f 

• r a| 


Rp is the feed location (3.75, 0, 0) 

Rr is the location of the center of the reflector (18.75, -32.5, 

-no) m an undeflected state. 

R A Is a unit vector in the direction of the reflector axis in 
Shuttle body coordinates 


The vector R A can be related to the direction cosine attitude matrices 
for the Shuttle body, Tj , and the reflector body, , by 



The relative alignment of the reflector to the Shuttle body is given by 
T 

T 1 T 4 which is a function of the structural deformations of the beam. 
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The line-of-sight error, e, is the angular difference between the 
target direction, given by the unit vector, D T , and the line-of-sight 
direction in Earth axes, T^qs* 

ARCSIN |D T X T iRl()S | or ARCSIN 

Computer programs are available which generate time histories of the 
rigid body and the mode shapes and frequencies for the body- beam-body 
configuration for "pitch" bending, "roll" bending and "yaw" twisting. 
Since the modes are based on solving explicitly the distributed parameter 
equations (without damping and without kinematic coupling) there is no 
limit to the number of modal characteristic sets that can be generated by 

the program. It will be the analyst’s decision as to how many modes need 
to be considered. 


Laboratory Experiment Description 

The second part of the design challenge is to validate in the 

laboratory, the system performance of the more promising control system 

designs of the first part. The experimental apparatus will consist of a 
dynamic model of the Space Shuttle orbiter with a large antenna reflector 
attached by means of a flexible beam. The dynamic model will be exten- 

sively instrumented and will have attached force and moment generating 
devices for control and for disturbance generation. A single, flexible 
tether will be used to suspend the dynamic model, allowing complete angular 
freedom in yaw, and limited freedom in pitch and roll. An inverted 

position will be used to let the reflector mast to hang so that gravity 
effects on mast bending will be minimized. The dynamics of the laboratory 
model will of necessity be different from the mathematical model discussed 
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earlier. 



Design Challenge, Part One 


For part one of the design challenge, the following mathematical 
problem is addressed. Given the dynamic equations of the Shuttle/antenna 

t t 

configuration, what control policy minimizes the time to slew to a target 
and to stabilize so that the line-of-sight (LOS) error is held, for a time, 
within a specified amount, 5, During the time that the LOS error is 
within 6, the attitude must change 90° to prepare for the next slew 
maneuver. This was previously referred to as the sescondary control task. 
The maximum moment and force generating capability will be limited. Advan- 
tage may be taken of selecting the most suitable initial alignment of the 
Shuttle/antenna about its assigned initial RF axis, line-of-sight. 

Random, broad band-pass disturbances will be applied to the configuration. 
Two proof-mass, force actuators may be positioned anywhere along the beam. 
The design guidelines are summarized below: 

1. The initial line-of-sight error is 20 degrees. 

e(o) = 20 degrees 

2. The initial target direction is straight down. 



3. The initial alignment about the line-of-sight is free to be chosen 
by the designer. Advantage may be taken of the low value of 
moment of inertia in roll. The Shuttle/antenna is at rest 
initially. 

4. The objective is to point the line-of-sight of the antenna and 
stabilize to within 0.02 degree of the target as quickly as 
possible. 

<5 =» 0.02 degree 

3 ?? 
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Cowro! moments can be applied .e ,00 Ha sampling r , t . eo ^ the 
shuttle and reflector bodies of 10,000 ft-lb for each axis. The 
commanded moment for each sals Is limited to 10,000 ft-lb. The 
actual control moment's response to the commanded value Is 
first-order with a time constant of 0.1 second. 

For the rolling moment applied to the shuttle body: 

-10 4 < R, , < io 4 

“* X, 1, command ~ 

M x,l (n + 1} = e " 0,1 \ 2 (n) + (1 - e -0 - 1 ) m (n) 

■ X ’ 1 X, 1 .command' 0 ' 

Equations for other axes and for the reflector body are similar. 
Control forces can be applied at the center of the reflector In 
the X and y directions only. The commanded force In a 
particular direction is limited to 800 lbs. The actual control 

force's response to the commanded value Is first-order with a 
response time of 0.1 second. 

For the side for applied to the reflector body: 


- 80 ° * F Y, command < «» 


F y (n + 1) = e -0 * 1 F (n) + (1 - e" 0 * 1 ) F (n) 

* Y, command'' 

Equations for X-axis are similar. 

Control forces using two proof-mass actuators (each having both 

X and y axes) can be applied at two points on the beam. The 

are limited to r 1 ft, and the masses weight 10 lbs each. 

The actual stroke follows a first-order response to limited 
commanded values* 
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For the X-axis of the proof-mass actuator at S 2 ^ 


-l < A y 0 . < 1 

— X, 2, command — 


A v 0 (n + 1) - e" 0 * 1 A (n) + (1 - e' 0 ' 1 ) A (n) 

X, 2, command 


X,2' 


Equations for other axes and locations are similar. 

8. The inertial attitude direciton cosine matrix for the Shuttle body 
lags in time the actual values by 0.01 second and are made at a 
rate of 100 samples per second. Each element of the direction 

cosine measurement matarix is contaminated by additive, 
uncorrelated Gaussian noise having an rras value of 0.001. The 
noise has zero mean. 



d n (n) 

d 12 (n) 

d 13< n) 

T s .measured^ 11 + ^ T s,true^ n ^ + 

d 21 (n) 

d 22 (n) 

d 23 ( n > 

where: 

d 31 (n) 

d 32 (n) 

d 33 (a) 

E{d ±j (n) } « 0 

E Kj (n)d kL (n) ^ = 0 

for i 

* k or 

j * L 

E Kj (n)d ij (n + k >i - ° 

for k 

* o 


= ( .001 J 2 

for k 

= 0 
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9. The angular velocity measurements for both the Shuttle and 
reflector bodies pass through a first-order filter with 0.05 sec 
time constant and lag in time the actual values by 0.01 second and 
are made at a rate of 100 samples per second. Each rate 
measurement is contaminated by additive, Gaussian, uncorrelated 

noise having an rms value of 0.02 degree per second. The noise 
has zero mean* 

For example: 


l,X>measured (n + ^ ,X, filtered (n) + e i,x (n) 

E K,X ( ") e i,x (n + k) } * 0 for k * 0 


- (.02)2 f or k = 0 

where 


1 ,X, filtered 


— 20 oo + 

1 ,X, filtered 


20 to 


1 ,X, true 


10. Three-axis accelerometers are located on the Shuttle body at the 
base of the mast and on the reflector body at its center. Two- 
axes (X and Y) accelerometers are located at intervals of 
10 feet along the mast. The acceleration measurements pass 
through a first-order filter with a 0.05 second time constant and 
lag in time the actual values by 0.01 second, and are made at a 
rate of 100 samples per second. Each measurement is contaminated 
by Gaussian additive, uncorrelated noise having an rms value of 
0.05 ft/sec^. 
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For example: 


a i ,X, measured^" + ^ a i ,X, filtered^ + T 1,X^ 

E ^ T 1 X (n) T 1 X (n + k) ) = 0 for k * 0 

=* ( .05)^ for k = 0 

where: 

a *i,X, filtered = " 20 3 1 ,X, filtered + 20 “l.X.true 

11. Gaussian, uncorrelated step-like disturbances are applied 
100 times per second to the Shuttle body in the form of 3-axes 
moments, having rms values of 100 ft-lbs. These disturbances 

have zero mean. 

For example: 

E ^ M D,X^ n) M D x < n + = 0 for k * 0 

= ( 100)^ for k = 0 

In summary, the designer's task for part one is to: (1) derive a 

control law for slewing and stabilization, coded in FORTRAN; (2) select an 
initial attitude in preparation for slewing 20 degrees; and (3) select two 
positions for the 2-axes proof-mass actuators. An official system 
performance assessment computer program will be used to establish the time 
required to slew and stabilize the Shuttle/antenna configuration. 
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Design Challenge, Part TVo 


As in part one, the task is to minimize the time to slew and stabilize 
a Shuttle/antenna configuration. The difference is that in part two of the 
design challenge, a physical laboratory model will be used instead of the 
dynamic equations of part one. The constraints on total moment and force 
generation capability will apply to part two, as for part one. Again, the 
analyst may select the initial alignment about the assigned initial RF 
line-of-sight. Disturbances will be injected into the Shuttle/antenna 
model. The designer's task will be similar to that for part one. 

CONCLUDING REMARKS 

A Design Challenge, in two parts, has been offered for the purpose of 
comparing directly different approach to controlling a flexible 
Shuttle/antenna configuration. The first part of the design challenge uses 
only mathematical equations of the vehicle dynamics; the second part uses a 
physical laboratory model of the same configuration. The Spacecraft 
Control Laboratory Experiment (SCOLE) program is being conducted under the 
cognizance of the Spacecraft Control Branch at the NASA Langley Research 
Center. The NASA/ IEEE Design Challenge has the advice and counsel of the 
IEEE-COLSS Subcommittee on Large Space Structures. Workshops will be held 
to enable investigators to compare results of their research. 
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MASS CHARACTERISTICS 











































Figure 2.- Schematic of the effect of bending on the 
1 ine-of-sight pointing error. 
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Roll bending mode ft 1. Frequency= .32 Hz 



Roll bending mode ft 2. Frequency= 1.29 Hz 



Roll bending mode ft 3. Frequency= 4.80 Hz 



Roll bending mode ft 4. Frequency= 12.29 Hz 



Roll bending mode ft 5. Frequency^ 23.68 Hz 



Roll bending mode ft 6. Frequency= 38.89 Hz 



Roll bending mode ft 1 . Frequency= 57.90 Hz 



Roll bending mode ft 8. Frequency= 80.72 Hz 
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Figure 4a.- Plots of normalized roll bending mode shapes 
for SCOLE configuration. 
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Pitch bending mode # 1. Frequency = .29 Hi 



Pitch bending mode # 2. Frequency^ 1.65 Hz 



Pitch bending mode # 3. Frequency= 497 Hz 



Pitch bending mode # 4 Frequency= 12.36 Hz 



Pitch bending mode # 5. Frequency= 23.72 Hz 



Pitch bending mode # 6. Frequency= 38.91 Hz 



Pitch bending mode # 7. Frequency^ 57.92 Hz 



Pitch bending mode # 8. Frequency= 80.73 Hz 


Figure 4b.- Plots of normalized pitch bending mode shapes 
for SCOLE configuration. 
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Torsional mode # 1. Frequency= .53 Hr 



Torsional mode # 2. 



Frequency= 45.12 Hz 
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Torsional mode # 3. Frequency^ 90.23 Hz 
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Torsional mode # 4. Frequency= 135.35 Hz 
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Torsional mode # 5. Frequency= 180.46 Hz 
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Torsional mode § 6. Frequency= 225.57 Hz 
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Torsional mode § 7. 



Frequency= 270.69 Hz 
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Torsional mode # 



8. Frequency= 315.80 Hz 
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Figure 4c.* Plots of normalized torsional mode shapes for 
SCOLE configuration. 


